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Abstract 

The Unitary Correlation Operator Method (UCOM) provides a means for nuclear structure 
calculations starting from realistic NN potentials. The dominant short-range central and 
tensor correlations are described explicitly by a unitary transformation. The application of 
UCOM in the context of the no-core shell model provides insight into the interplay be- 
tween dominant short-range and residual long-range correlations in the nuclear many -body 
problem. The use of the correlated interaction within Hartree-Fock, many-body perturba- 
tion theory, and Random Phase Approximation gives access to various nuclear structure 
observables throughout the nuclear chart. 
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1 Introduction 



In recent years several realistic nucleon-nucleon interactions like the Argonne VI 8 
[1] and the CD Bonn potentials as well as interactions derived from a chiral effec- 
tive field theory [2,3] have been constructed on the basis of high-precision nucleon- 
nucleon scattering data. These potentials are used in ab initio nuclear structure cal- 
culations throughout the p-shell, e.g., in the framework of the Green's function 
Monte Carlo method or the no-core shell model [4]. The use of these realistic po- 
tentials for nuclear structure studies in heavier nuclei poses an enormous challenge. 
Traditional many-body methods, like the Hartree-Fock approach or the Random 
Phase Approximation, cannot be used in connection with the bare NN interaction. 
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The reason is the inability of the simple model spaces underlying these approaches 
to describe the dominant short-range correlations, which are present in the exact 
many-body eigenstates. 

The two most important types of many-body correlations are those induced by the 
short-range repulsion and the strong tensor part of the NN interaction. Already in 
the deuteron their presence is evident: (/) The probability amplitude for finding the 
two nucleons at very small distances is strongly depleted as a result of the short- 
range repulsive interaction, (ii) Apart from the L = component the ground state 
exhibits a L = 2 admixture which is essential for the binding and is generated by 
the tensor part of the NN interaction. Also for heavier nuclei, these central and ten- 
sor correlations have a dominant impact on the structure of the many-body state. 
Neither of these correlations can be described properly by a single or a superpo- 
sition of few Slater determinants. Therefore, a naive inclusion of the bare realistic 
NN-potential into a Hartree-Fock-type calculation has to fail. 

For nuclei beyond the p-shell one is bound to use simplified model spaces for an 
approximate solution of the many-body problem. Therefore, the short-range cor- 
relations have to be accounted for explicitly, e.g., by transforming the bare real- 
istic interaction into an effective interaction adapted to the available model space. 
One possible approach to construct a phase-shift equivalent effective interaction 
is the Unitary Correlation Operator Method discussed in Section 2. The result- 
ing correlated interaction V UC om can then be used as a universal input for differ- 
ent many-body approaches, ranging from the no-core shell model (Section 3) over 
Hartree-Fock and many-body perturbation theory (Section 4) to the Random Phase 
Approximation (Section 5). 



2 Unitary Correlation Operator Method 



The basic idea of the Unitary Correlation Operator Method (UCOM) is to include 
the dominant correlations into the many-body state by means of a unitary transfor- 
mation [5,6,7,8]. Starting from an uncorrelated many -body state in the sim- 
plest case just a Slater determinant, a correlated state |\I>) is defined through the 
application of the unitary correlation operator C: 

|§) = C|W). (1) 

Alternatively, one can perform a similarity transformation of the operators of all 
relevant observables (e.g. the Hamiltonian, coordinate and momentum space den- 
sities, transition operators, etc.): 

6 = C f OC . (2) 



2 



Due to unitarity both approaches are equivalent. For most many-body calculations 
the formulation through correlated operators is, however, more convenient. 



The correlation operator C is decomposed into a central correlator C r and a ten- 
sor correlator Cn, reflecting the two dominant types of correlations in the nuclear 
many-body problem: 



C = CnC r = exp 



i<3 



iSgn(u) exp -i^g r (ij) 



i<3 



(3) 



Both operators are defined as exponentials of Hermitian two-body generators 
and g r , respectively. They are given in a closed analytic form which reflects the 
mechanism by which correlations are induced by the interaction. 

The task of the central correlator C r is to generate the hole in the two-body density 
distribution at small particle distances caused by the repulsive core in the central 
part of the interaction. Pictorially speaking, C r has to shift pairs of particles that 
are closer than the core radius apart from one another. The two-body generator 
for this distance-dependent shift can be written as g r = |[s(r)q r + q r s(r)], where 
q r = |[q • (r/r) + (r/r) • q] is the radial component of the relative momentum q of 
a particle pair. The function s(r) determines the distance-dependence of the shift. 
It is large for small r and vanishes at large distances. 

The tensor correlation operator Cq has to generate the complex entanglement be- 
tween the angular structure of the relative two-body states and the spin orientation. 
An essential ingredient is the component of the relative momentum q perpendicular 
to r, the so-called orbital momentum q n = q — ^ q r . The generator has the form 

En — f$(r) (<x r q n ) (er 2 - r) + (r <-> q n ) which is similar to the tensor operator Si 2 . 
The function $(r) describes the magnitude of the shift as a function of distance. 

For the following many-body calculations, the notion of correlated operators is ad- 
vantageous. The operators of all observables under consideration have to be trans- 
formed consistently. Since the correlation operators are defined as exponentials of 
two-body operators, the correlated operators contain irreducible contributions for 
all particle numbers. We organize the different irreducible terms according to their 
rank in a cluster expansion 

H = C f HC = H [1] + H [2] + H [3] + • • • . (4) 



Here we used the Hamiltonian H = T + V as an example, but the same holds true 
for any other operator. If the range of the correlators is sufficiently small compared 
to the average particle distance in the many-body system, three-body and higher 
order terms in the cluster expansion are small and we can restrict ourselves to the 
two-body approximation 

gC2 = + f[2] + y[2] = T + VucQM ; (5) 
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where T M = T and T' 2 ' are the one- and two-body contributions of the correlated 
kinetic energy, resp., and V'^ is the two-body part of the correlated NN-potential. 
All two-body contributions are subsumed in the correlated interaction Vucom- It is 
by construction phase-shift equivalent to the original, uncorrected NN-potential as 
long as the correlators have finite range. 

The remaining task is the determinantion of the correlation functions s(r) and 
i?(r) entering into the generators of the unitary transformations. For each spin- 
isospin channel their parameters can be obtained from an energy minimization in 
the two-body system. This procedure and the optimal correlators for the Argonne 
VI 8 (AVI 8) potential are discussed in Ref. [8]. The tensor correlation functions re- 
quire a special treatment. Since it originates from the one-pion exchange, the tensor 
force is long-ranged, and so are the tensor correlations induced in the two-body sys- 
tem. In a many-body system, the long-range component of the tensor correlations 
between two nucleons is screened due to the presence of other nucleons. In antici- 
pation of this effect, we restrict the range of the tensor correlator by a constraint on 
the integral of the correlation function, 1$ — J drr 2 i!)(r). Hence, only short-range 
correlations are described explicitly by the unitary transformation. Long-range cor- 
relations have to be covered by the many-body states — this will be illustrated in the 
following sections. 



3 No-Core Shell Model Calculations 

As a first application of the correlated realistic interaction V UC om we consider a 
straightforward no-core shell model diagonalization within a harmonic oscillator 
basis. The shell model basis itself is able to describe part of the many-body cor- 
relations, depending on the size of the model space. Hence the dependence of the 
energy on the model-space size provides information on the role of short-range 
correlations and on the contribution from residual long-range correlations. For the 
calculations we employ the translationally invariant no-core shell model code de- 
veloped by Petr Navratil [9], but without using the Lee-Suzuki transformation. The 
computation of the relevant two-body matrix elements of Vucom in the harmonic 
oscillator basis and further results are discussed in Ref. [8]. 

Figure 1 shows the ground state energy of 4 He as a function of the oscillator pa- 
rameter Ml, for different sizes of the model space, characterized by the maximum 
relative oscillator quantum number jV max . The upper panel corresponds to a calcu- 
lation with the bare AVI 8 potential. Evidently, even for the largest feasible model 
spaces, the energy is not yet converged. The reason is that a full description of short- 
range central and tensor correlations requires even larger model spaces, which are 
computationally not tractable. The picture changes if we use Vucom, include 
the unitary transformation of the Hamiltonian. The convergence is dramatically 
improved since the short-range central and tensor correlations are now treated ex- 
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Fig. 1. Results of no-core shell model calculations using the correlated AVI 8 potential. Left 
panel: convergence of the ground state energy of 4 He for bare (upper plot) and correlated 
AV18 potential (lower plot). Right panel: Tjon-line and dependence of the energy on the 
correlator range as described in the text (taken from [8]). 

plicitly by the unitary correlation operator. Note that a bound nucleus is already 
obtained with a single Slater determinant (i.e. W max = 0). With increasing size of 
the model space, the ground state energy is lowered further. This is the result of the 
improved description of long-range correlations — not accounted for by the unitary 
transformation — by the model space. 

A second interesting aspect is illustrated on the right-hand side of Fig. 1, where 
the converged ground state energies of 3 H and 4 He are plotted. Each data point 
corresponds to a different interaction. The exact energies for the different bare NN- 
interactions, like the AVI 8, the CD Bonn and the Nijmegen interactions (circles), 
fall onto the so-called Tjon-line [10] but are far away from the experimental point. 
Three-nucleon interactions (diamonds) are needed to obtain binding energies in 
the experimental region. The exact energies for the correlated interaction Vucom 
based on AVI 8 (triangles) depend on the range 1$ of the triplet-even tensor corre- 
lation function. With increasing range the energy is lowered and the full Tjon-line 
is mapped out. This is related to the omission of three-body (and higher-order) 
terms in the cluster expansion of the correlated Hamiltonian. If these terms were 
included, the energies would be exactly the same, independent of the correlator 
range, because of the unitarity of the transformation. The fact that the range of the 
tensor correlator can be chosen such that the energies are close to experiment (e.g. 
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Fig. 2. Ground state energy of various closed shell nuclei obtained with the correlated AVI 8 
potential within a HF calculation (circles) and in HF + MBPT (squares and diamonds) in 
comparison to experimental binding energies (bars) (taken from [11]). 

for h = 0.09 fm 3 ) can be explained by a cancellation between genuine three-body 
forces and the induced three-body contributions of the cluster expansion. In other 
words, the impact of the net three-body force on the binding energies can be mini- 
mized by a proper choice of the correlator range. 



4 Hartree-Fock & Many-Body Perturbation Theory 



Using the Vucom interaction fixed within the no-core shell model we perform 
Hartree-Fock (HF) calculations of nuclear ground states throughout the nuclear 
chart. Since the HF many-body state (Slater determinant) alone is not able to de- 
scribe any many-body correlations, the use of bare realistic interactions does not 
lead to bound nuclei. The explicit inclusion of the short-range correlations, e.g., via 
the unitary correlation operators is inevitable. 

We have implemented the HF scheme in the harmonic-oscillator representation, 
using the translationally invariant Hamiltonian H int = T — T cm + Vucom* where 
Vucom contains charge-dependent and Coulomb terms [11]. The results for ground 
state energies of closed-shell nuclei ranging from 4 He to 208 Pb are depicted in Fig. 
2. The optimal correlator for 1$ = 0.09 fm 3 is used, and the single-particle basis 
includes 13 major oscillator shells. Evidently, the HF binding energies are signif- 
icantly smaller than the experimental ones. This is not surprising, since residual 
long-range correlations as they appeared in the no-core shell model calculations 
cannot be described by the HF ground state. 

An estimate for the impact of residual long-range correlations on the binding ener- 
gies can be obtained within many-body perturbation theory. The evaluation of the 
second and third order perturbative contributions on top of the HF result is straight- 
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forward [11]. Figure 2 summarizes the results for the ground state energies includ- 
ing second order correlations (for light nuclei also third order). Again, 13 major 
oscillator shells are included to obtain a satisfactory degree of convergence for the 
perturbative correction. The agreement with the experimental binding energies per 
nucleon is remarkably good throughout the whole mass range. The absence of any 
systematic deviation for larger mass numbers proves that the cancellation between 
genuine three-body force and induced three-body contributions, which we observed 
in the no-core shell model for light isotopes, works throughout the nuclear chart. 
Furthermore, the calculations establish the perturbative character of the long-range 
correlations. Note that a perturbative treatment of the short-range correlations is 
not possible — in our approach they are covered by the unitary correlation operators 
from the outset. 

However, the good agreement with experimental data does not hold for all observ- 
ables. The charge radii obtained in HF for heavier nuclei are too small in com- 
parison to experiment [11]. The inclusion of perturbative corrections improves the 
result but still leaves deviations of up to 1 fm for the heaviest nuclei. This is an 
indication that a net three-body force is needed to reproduce all observables, al- 
though its impact on the energy might be small. This issue is the topic of future 
investigations. 



5 Random Phase Approximation 

In addition to global ground state properties, collective excitations provide a value- 
able probe to understand the role of correlations in the nuclear many-body problem. 
We use the standard Random Phase Approximation (RPA) [12] and its extensions 
[13] to study the behavior of collective excitations based on V UC om- Starting from 
the HF solution for the ground state we solve the RPA equations in a fully self- 
consistent way using the same intrinsic Hamiltonian as for the HF treatment. In 
this way the spurious center-of-mass mode appears fully decoupled at very low 
excitation energies of the order of lOkeV and the energy-weighted sum rules are 
fulfilled with maximum deviations of ±3%. 

The results obtained for isoscalar monopole, isovector dipole, and isoscalar quadru- 
pole excitations in 40 Ca, 90 Zr, and 208 Pb using the standard Vucom interaction are 
summarized in Figure 3. In all cases a collective resonance appears in the RPA re- 
sponse, which is not trivial since we use a realistic NN-interaction. The centroid 
energies for the isoscalar giant monopole resonances are in nice agreement with 
experiment. Keeping in mind that there are no free parameters, this is a remark- 
able result indicating that the incompressibility generated by Vucom is reasonable. 
However, for the isovector dipole and the isoscalar quadrupole giant resonances the 
calculated centroid energies are systematically larger than the experimental ones. 
This hints at too small a value for the effective mass which is consistent with the 
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Fig. 3. RPA strengths distributions for isoscalar monopole (ISM), isovector dipole (IVD), 
and isoscalar quadrupole (ISQ) excitations in 40 Ca, 90 Zr, and 208 Pb using Vucom- The 
curves result from a Lorentzian folding of the discrete strength distributions and the black 
triangles indicate the centroid energies extracted from experiment (see also [12]). 

very wide single-particle spectra resulting from the HF calculations [11]. Again 
this could be an indication of a residual repulsive three-body force but also of the 
importance of additional correlations not included in the standard RPA framework. 



6 Conclusions & Outlook 



We have used the correlated realistic interaction Vucom derived from the Argonne 
VI 8 potential as a universal starting point for nuclear structure calculations through- 
out the nuclear chart. Different many-body approaches ranging from no-core shell 
model to Hartree-Fock, many-body perturbation theory, and RPA have been em- 
ployed using the same interaction. We observe that the binding energies per nu- 
cleon resulting from Vucom are in good agreement with experiment through the 
whole mass range, indicating that the impact of residual three-body forces on this 
observable is minimal. However, for other observables systematic deviations from 
experiment emerge minly in heavier nuclei: (i) the rms-radii are too small, (ii) the 
mean level-spacing of the HF single-particle spectra is too large leading to a too 
small effective mass, {Hi) the centroid energies of isovector dipole and isoscalar 
giant quadrupole resonances are overestimated. One possible origin of these devia- 
tions, besides missing long-range correlations, is a residual repulsive three-body in- 
teraction not included in the present calculations. Initial calculations using a simple 
phenomenological three-body contact interaction indicate that all aforementioned 
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discrepancies can be reduced. The inclusion of effective three-body interactions 
into the many-body schemes discussed here is the subject of ongoing investiga- 
tions. 
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